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ABSTRACT‏ 
A STUDY OF BOILING WATER FLOW REGIMES AT LOW PRESSURE‏ 


Mario P. Fiori 


Submitted to the Department of Nuclear Engineering in partial 
fulfillment of the requirement for the joint degrees of 
Nuclear Engineer and Master of Science in Mechanical Engineer- 
Ange 


A comprehensive experimental program to examine flow 
regimes at pressures below 100 psia for boiling of water in 
tubes was carried out. 


An electrical probe, which measures the resistance of 
the fluid between the centerline of the flow and the tube 
wall, was used to identify the various flow regimes. This 
probe proved to be an ideal detection device, because of its 
simplicity, reproducibility, and accurate representation of 
the flow pattern within the heated test section. 


The major flow regimes observed were bubbly, slug and 
annular flow. Under certain conditions at high flow rates, 
a wispy-annular flow pattern was observed. The effects of 
mass velocity (0.2 x 100 - 2.4 x 100 1bm/nr-1t2), inlet temp- 
erature 0 150, 200°F), exit pressure (30, 100 psia), 
quality (x = -10 - +7 percent), purity (9, 40 PPM NaCl; 1-3 
megohm-cm), length (L/D=30, 60, 90), diameter TT 0.242 
in.), and orientation (vertical and horizontal) on the flow 
regimes were studied. Flow regime maps on coordinates of 
mass velocity and quality are presented for these conditions. 


Bubbly and slug flow occurred primarily in the sub- 
cooled region, while fully developed annular flow was reached 
at equilibrium qualities between 2 and 4 percent. The tran- 
sitions between the different flows were shifted to regions 
of increased subcooling when velocity, pressure, and heat 
flux increased, and when inlet temperature decreased. Purity 
and geometry had little affect on the flow regime boundaries. 


The shifting of the transitions is related to the agglom- 
eration point, which is that point at which the bubbles so 
coalesce that slug flow is first observed. The agglomeration 
point depends on the point of incipient boiling, the number 
of bubbles in the flow, and the number of collisions per 
bubble. These latter quantities in turn depend on velocity, 
temperature, pressure, and heat flux. 


The flow regime information obtained in this study should 
be of value in correlating and interpreting low pressure heat- 
transfer data. The flow regime data were found to be useful 
in explaining the effect of inlet temperature on burnout heat 
1 , 
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NOMENC LATURE 
heat transfer area (ft?) 
specific heat (Btu/lbm°F) 
diameter (in.) 
mass velocity (1bm/hr-ft*) 
heat transfer coefficient (Btu/hr-ft<oF) 
enthalpy at exit (Btu/lbm) 
enthalpy of saturated liquid (Btu/lbm) 
latent heat (Btu/lbm) 
enthalpy at entrance to test section (Btu/lbm) 
shunt voltage (mv) 
thermal conductivity (Btu/hr-ft°F) 
length (in., ft) 
pressure (psia) 
pressure at exit of test section (psia) 
rate of heat flow (Btu/hr) 


heat flux (Btu/hr-ft^) 


= heat flux at incipient boiling (Btu/nr-1tó) 


critical heat flux (Btu/hr-ft^) 
oscilloscope internal impedance (ohms) 
resistance between probe and tube wall (ohms) 
variable resistance (ohms) 

water temperature at inlet or exit of tube (°F) 
water bulk temperature (°F) 

wall temperature (°F) 


saturation temperature (°F) 


voltage across test section (volts) 
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Vz = battery voltage (volts) 

V+ = voltage across variable resistor (volts) 
w = mass flow rate (lbm/hr) 

x = quality 

u = viscosity (lbm/hnr-ft) 


Flow Regime Designations 


B = @ = bubbly 
BTS = & = bubbly-to-slug transition 
S = Q = slug 
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STA = slug-to-annular transition 
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Chapter I 
INTRODUCTION 
1.1 Background of Problem 

In recent years the interest in two-phase, gas-liquid 
flow has dramatically expanded as is indicated by the tre- 
mendous increase in the number of experiments and published 
papers in this field. If one examines the literature, one 
will find that a relatively large number of models have 
been proposed for analyzing and correlating experimental 
pressure-drop, heat-transfer, critical heat flux, and void- 
fraction data. These models frequently depend on the flow 
regimes (or flow patterns) which are geometric descriptions 
of the distribution of the liquid and gas phases within a 
channel. qouse(1)* considers knowledge of the flow regime 
to be as important as knowing whether the flow is laminar 
or turbulent in single-phase flow. 

Baker (2) showed that it was not practical to disregard 
the flow regimes in devising a two-phase pressure-drop 
correlation. Martíne111 1513) assumptions in his pressure- 
drop correlations tend to limit them to annular flow. 

Lavin and Young propose boiling heat-transfer correla- 
tions for refrigerants, which are dependent on the flow 
regimes. Numerous models for the critical heat flux condi- 
tion are based on the postulate of spray annular flow; 
however, there are some analyses which assume that fog flow 


is present. 


“Numbers in parenthesis refer to References on p. 80. 
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It 1s apparent, then, that it is important to know 
the flow regime 1f one is to utilize an appropriate model; 
however, direct relationships between flow regimes and 
Obsérvables, i.e., pressure drop, heat-transfer coeffi- 
Client, critical heat flux, and void fraction have not been 
clearly established. U^? Itigulll appears, though, that 
flow regime observations will be helpful in reducing the 
complications of correlating the quantities of basic 
interest. 

A survey of the literature indicates that flow regime 
studies have almost always been conducted under adiabatic 
conditions. However, exploratory studies have shown that 
flow patterns are considerably altered when there is heat 
addition, particularly in low pressure tens (© 
Comprehensive flow regime maps have been prepared at 
Dynatecn(7)(8) and Harpe lil”) for high pressure boiling 
water systems. No similar study has been conducted at low 
pressures (below 100 psia). 

Low pressure systems are receiving increased atten- 
Cion, particularly in the nuclear reactor and desaliniza- 
tion programs. Ina recent study of critical heat fluxes 
at low pressure, Lopina (2°) observed a pronounced effect 
of inlet temperature. For constant geometry, flow rate, 
and pressure level, the critical heat flux increased as 
the inlet temperature increased. This effect, which is 


contrary to what is observed at high pressures, appears to 


be related to the flow regimes. The phenomenon of 


c 
upstream boiling burnout in uniformly heated channels also 
appears to be related to flow regimes, in particular, the 


slug flow regime. 


1.2 Flow Regimes 


There have been various descriptions of flow regimes 
by different investigators. These classifications have 
often been quite arbitrary, and are also complicated by the 
fact that vertical and horizontal flows cause different 
regimes. Baker (2) described no less than seven different 
regimes in adiabatic two-phase, gas-liquid flow in horizon- 


(11) and Kepple and Tung (22) have presented 


tal pipes. Vohr 
complete bibliographies on studies of adiabatic flow 
patterns. 

The flow regimes which may be encountered in flows in 
tubes with vertical upflow can be described in the order of 
increasing gas or vapor content. If very little gas or 
vapor is involved, the flow will be characterized by isola- 
ted vapor or gas bubbles dispersed rather uniformly in the 
continuous liquid phase. This is known as "bubbly flow". 
More bubbles appear with an increase of vapor or gas until 
the bubbles begin to touch one another and agglomerate. 

As more bubbles agglomerate, the gas begins to flow in 
long cylindrical, bullet-shaped bubbles, Between the 
bubbles are slugs of liquid and hence "slug flow", 

Further increase of the gas content destroys the stability 
of the bubbles, producing "churn flow" or "semi-annular 


flow" in which slugs of liquid separate adjacent gas 
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bubbles only temporarily. As the gas content continues to 
increase, the liquid flows on the wall as an annular film 
and the gas or vapor flows in the center of the tube. This 
geometry is called "annular flow". "Fog flow" or "dispersed 
flow" then follows where the gas or vapor phase occupies 
the entire tube but small drops of liquid are dispersed 
within it. This flow is also called "mist flow". 

Other flow patterns which have been encountered are 
"froth flow", “spray annular flow" and "wispy annular flow". 
Froth flow is quite similar to bubbly flow except the 
bubbles are so numerous that the flow has a frothy appear- 
ance. Spray-annular flow is a combination of mist and annu- 
lar flow. Some of the liquid flows in an annular film and 
some flows as small droplets in the vapor core. Wispy annu- 
lar is similar to spray annular except the flow in the core 
consists of agglomerated liquid instead of small droplets. 

There are very similar flow descriptions for horizontal 
flow. However, stratification is sometimes encountered in 
horizontal flow and several new regimes are introduced. 
"Plug flow" is similar to slug flow except that the bubbles 
move along the upper part of the channel or pipe. In 
"stratified flow" the liquid flows at the bottom, and gas or 
vapor at the top of the pipe. 

These regimes appear to be quite well defined. However, 
Lor practical reasons it is desirable to minimize the 
number of regimes. Gouse (+) suggests that there are only 


four distinct flow patterns for vertical upflow: bubbly, 
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slug, annular, and mist; while all the others which have 
been observed are merely transitions from one to another. 

The above flows were described for adiabatic condi- 
tions. In heated systems such clean geometric descriptions 
of the flow are not expected because of several character- 
istics not encountered under adiabatic conditions. Some of 
these characteristics are: 

1) Vapor is continually formed along the tube. This 
creates a volume flow rate which continuously increases 
along the tube; thus the flow regimes would be expected to 
be different in different sections of thetube. 

2) Bubbles may be generated and disturb the flow 
pattern. 

In any case, the adiabatic flow regimes are a basis for 
interpreting and classifying the flow regimes with heat 
addition. The classification and both qualitative and semi- 
quantitative description of the flow patterns observed in 


this study are discussed in Section 2.3.3. 


1.3 Flow Regime Detection 
Numerous methods of studying flow regimes in heated 


systems have been developed. High-speed still and movie 
photography has been most popular. More exotic methods 
which have been used include isokinetic flow sampling, 
gamma or beta ray attenuation, and hot-wire anemometry. 
Each of these methods has drawbacks. Either there is too 
much data to analyze, as in the case of movies, or the 


methods are very complicated and expensive. The electri- 





Eos 


cal resistance probe which was developed by Solomon i3) 
specifically for flow regime detection, is an inexpensive, 
simple, and reliable device. 

This probe measures the resistance between the center- 
line of the fluid flow and tne tube wall. A bridge of 
liquid across the tube will cause a closed circuit, and 
gas or steam across the tube will result in an open circuit. 


(14) and Haberstroh(15) have successfully used this 


Nassos 
type of probe for adiabatic flow. Griffith "® adapted 
this probe for a high pressure boiling water system and 
delineated the slug-to-annular transition. Dynatecn / (9) 
has recently been using this technique to map boiling water 


flow regimes at high pressures (500-1000 psia). 


1.4 Scope of Research 


A comprehensive study to examine boiling water flow 
regimes in tubes of small diameter at pressures below 100 
psia was undertaken. A resistance probe described in 
detail in Section 2.2 was chosen as the primary means of 
flow regime detection; however, a visual section was also 
employed. The experimental program consisted of an inves- 
tigation of the effects of velocity, inlet temperature, 
exit pressure, quality, (heat flux), purity, length, diam- 
eter, and orientation on the flow regimes. The first 
phase of the study was devoted to probe development and 


interpretation. 
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Chapter II 


EXPERIMENTAL PROGRAM 


21 Description of Apparatus 
21 上 Hydraulic System 


The experimental facility used was the low-pressure 
test loop located in the MIT Heat Transfer Laboratory. The 
basic apparatus was designed and constructed in 1961.07) 

A schematic of the loop is presented in Fig. 1. The 
pipings and fittings are made of brass and stainless steel 
for corrosion resistance. Rayon reinforced rubber hose is 
used where flexible connections are required. Distilled 
water is circulated by a bronze, two-stage, regenerative 
pump providing a discharge pressure of 260 psig at 3.6 gpm. 
The pump is driven through a flexible coupling by a 3-hp 
Allis-Cnalmers induction motor. A Fulflo filter is in- 
stalled at the pump inlet. Pressure fluctuations at the 
outlet of the pump are damped out by means of a 2.5-gal 
Greer accumulator charged with nitrogen to an initial 
pressure of 40 psig. This accumulator contains a flexible 
bladder-type separator which prevents the nitrogen from 
being absorbed by the system water. After the accumulator, 
the flow splits into the by-pass line and the test-section 
lines 

In the test-section line, fluid flows through a 
Fischer-Porter flowrator followed by a preheater, thence 
through a Hoke metering valve and the test section, after 
which it merges with fluid from the by-pass line. The 





4! 
= 





-Ba 
flow then goes through the heat exchanger and returns to 
the pump. The preheater consists of four Chromalox heaters 
of approximately 6 kw each. Three of these are controlled 
simply with "off-on" switches while the fourth can provide 
a continuous range from O to 6 kw by means of a bank of 
two variacs mounted on the test bench. Quick-action James- 
bury ball valves are installed before the inlet to the 
flowrator and after the exit from the test section. This 
permits quick isolation of the test section to minimize 
fluid loss when conducting burnout tests. The exit valve 
is also used to adjust the test-section pressure. 

Flow through the by-pass line is controlled by a ball 
valve on each side of which there is a 300-psig pressure 
page, Pump operating pressure, and hence the pressure up- 
stream of the test section, is controlled by this valve. 

The heat exchanger is a counterflow type with system 
water flowing in the inner tube and city water in the 
outer annulus. Except at very high power levels, the heat 
exchanger maintained a constant pump inlet temperature. A 
Fulflo filter is installed on the city water line to 
reduce scale formation in the exchanger, 

The distilled water was de-ionized continuously, 
except during contamination studies, by passing a portion 
of the flow through four mixed-bed resin demineralizer 
cartridges installed in a Barnstead- "Bantam" demineralizer 
unit. A 4.7-gal degassing tank is provided with five elec- 


trical heaters (3-200 vac and 2-110 vac). This tank also 
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serves as a surge tank. A 15-gal stainless-steel storage 
tank for filling the system is mounted directly above the 
degassing tank and can be filled with distilled water from 


standard 5-gal bottles with a small Hypro pump. 


2.1.2 Power Supply 


Test-section electrical power is provided by two 
36-ku de generators connected in series. Each generator is 
rated at 12 volts and 3000 amperes. The power control con- 
sole permits coarse or fine control from 0 to 24 volts. A 
water-cooled shunt installed in parallel with the test sec- 
tion protects the generators against the shock of the sudden 
open circuit which occurs at burnout. Power is transmitted 
from the main bus to the test section by water-cooled power 
leads. In the present tests, the downstream or exit power 
lead was at ground potential. Rubber hose connected both 
inlet and exit chamber plenums to the main loop in order 
to electrically isolate the test section. The exit of the 
test section was then separately grounded (see schematic 


of test section, Fig. 2). 


2.1.3 Instrumentation 

All temperatures were measured by copper~-constantan 
thermocouples made from 30-gauge Leeds and Northrup duplex 
wire. The test-section inlet temperature was measured by 
a thermocouple directly in the fluid stream, upstream of 
the Hoke metering valve. The thermocouple was introduced 
at this location through a Conax fitting equipped with a 


lava sealant. 
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All pressures were read on Bourdon-type gages located 
as shown in Fig. 1. The test-section inlet and exit pressures 
were measured with Helicoid 8-1/2 in. gages of 200 psig and 
100 psig, respectively. Both are specified to an accuracy 
of +0.25% of full scale. 

A are of metering tubes and floats, which could be 
installed interchangeably in the basic Fiscner-Porter flow- 
meter housing, provided measurement of the test-section flow 
from 1.5 to 4000 lbm/hr. 

The voltage drop across the test section was read 
directly on a Weston multiple-range dc voltmeter with a spec- 
ified accuracy of 4 1/2 percent. The current flow was 
determined by using the Minneapolis-Honeywell, Brown recor- 
der to measure the voltage drop across a calibrated shunt 
(60.17 amp/mv) in series with the test section. 

The details of the resistance probe are discussed in 
section 2.2. Polaroid cameras were used for the photo- 
graphic study. The camera settings and procedures are dis- 


cussed in Appendix B. 


2.1.4 Operating Procedures 
After the test section was installed, the loop and de- 


gassing tank were filled with distilled water from the 
supply tank. The water in the degassing tank was then 
brought to a boil while the loop water was circulated with 
the heat-exchanger coolant off. Degassing was accomplished 
by by-passing a portion of the loop water into the top of 


the vigorously boiling degassing tank. This was continued 
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until the temperature of the loop rose to approximately 
180°F. A standard Winkler analysis described in Ref. 17 
indicated that this method of degassing reduced the air 
content to less than 0.1 cc air/liter. After completing 
degassing, the loop water was adjusted. to the desired in- 
let temperature and the desired pressures were set. A 
pressure of from 210 psig to 260 psig was maintained at all 
times on the upstream side of the Hoke metering valve to 
prevent system-induced instabilities. The generators were 
then started and allowed to warm up. 

Nine burnout tests were made in order to have reliable 
critical heat flux data for conditions similar to the planned 
flow regime studies. This knowledge aided in avoiding burn- 
out of the probe test sections. Power was applied to the test 
section in small increments while maintaining a constant flow 
rate and inlet temperature. The test-section flow rate, 
shunt voltage (test-section current), inlet temperature, and 
exit pressure were recorded for each voltage setting. The 
critical or burnout heat flux was defined to be that heat flux 
which caused physical failure of the test section. After 
commencing flow regime studies, it was seen that the area of 
interest was usually far removed from (A/a) one Consequently 
no further burnout studies were made. Appendix A contains a 
tabulation of the burnout points. 

The same data were recorded for the flow regime runs. 

A run consisted of taking data at a particular mass flow 


rate for a given set of conditions, i.e. D, L/D, Ty» Pe» 
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As the power was increased, various responses on the oscil- 
loscope indicated the type of flow (see Probe Interpretation, 
Section 2.3). Water purity was checked by means of a 
Barnstead Purity Meter (0.1 - 18 megohm-cm) or an Industrial 


Instruments Solu-Bridge Meter (O - 30 PPM NaCl). 


2.2 Resistance Probe and Test Section 

2.2.1 Principle of Operation 

The resistance probe was used to determine the flow 
regime at the exit of the test section. The probe is a wire 
placed in the center of the tube as shown in Fig. 2. This 
wire is covered by teflon sleeving except for the tip of the 
wire. The circuit diagram for the probe is shown in Fig. 3a. 
When there is a bridge of water across the tube at the test- 
section exit, there is essentially a closed circuit between 
the probe and the tube wall. This results in a maximum volt- 
age across resistor Ry which is measured by an oscillo- 
scope. On the other hand, steam passing between the probe 
and the wall resembles an open circuit giving a voltage 


reading across Ry near zero. 


2.2.2 Probe Integration with Test Section 
Two types of test sections were used for this study. 


The first type was used for preliminary burnout studies. 
The second type incorporated the resistance probe and a 
sight section. The sight section was not always used 
because the rubber gasket would not hold for exit pressures 


above HO psig. The heated sections were made of standard 
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304 stainless tubes. The burnout test sections and the first 
probe test. section had an 0O.242-in. 1.d. and 0.3125=1n. 0.0. 
while the second probe test section had an i.d. of 0.094 in. 
and o.d. of 0.120 in. The details of the two test sections 
are shown in Fig. 2. 

In the burnout test section, 3/4-in. brass bushings 
were silver-soldered to the stainless steel tube to form 
the power connections upstream and downstream. A pressure 
tap was installed in each exit bushing with a No. 80, 
0.0135-in. diameter hole drilled through the tube wall. 

The pressure tap was placed here rather than in the exit 
plenum chamber because of the choking phenomena whereby the 
plenum pressure is lower than. the exit pressure , (10) 

The probe test section had pressure taps located at 
both the inlet and exit. The exit section was modified to 
include the probe and pyrex sight section. The inside 
diameter of the glass tubing was equal to the outside diam- 
eter of the test section. This caused a slight expansion 
at the entrance to the visual section. 

The brass supporting structure for the probe, pressure 
tap, and sight section also served as the exit power connec- 
tion. The pressure tap is similar to the previous one 
described. 

The probe was fabricated from 0.035-in. diameter temp- 
ered stainless-steel spring wire for the 0.242-in. diameter 
tube and from 0.015-in. wire for the 0.094-in. diameter 


tube. Teflon spaghetti was slipped over this wire and 


۷ 
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stretched until the teflon set snugly on the wire. The 
final probe diameters were 0.060 in. and 0.031 in. Even the 
small diameter probe had adequate stiffness to prevent vi- 
bration which could create spurious signals. 

The straight probe was pushed into the test section 
from the packing gland side until it touched the opposite 
wall. The tip was then carefully bent and the teflon was 
forced around the wire until only a 1/4 in. of probe tip was 
bare. 

The teflon acts as an insulation between the probe and 
the tube wall. More important, since teflon is not wetted 
by water, there is no chance of a short-circuit between the 
exposed end of the wire and the tube wall. The drying time 
of the teflon after passage of a slug of liquid is extremely 


9 


short, less than 10 " seconds as observed with an oscillo- 


agape. t2) 


2,2,3 Probe Sensitivity 


A schematic diagram of the probe circuit is shown in 
Fig. 3a. It can be shown that the voltage reading, Vi? is 
related to the resistance of the water-steam mixture in the 


test section, Ry by the equation 


- Yg (1) 


í 1 ۸ m 7 — 
Rul” + RO) + 1 


1 


Furthermore, it can be shown that the maximum sensitivity, 


aV 
defined by a maximum in (SE) is obtained when 
۷٧ 


1 1 1 
a sa pa (2) 
Ru Ry Ry 
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At this maximum sensitivity, the instrument voltage, V is 


Ie 
related to the battery voltage, Vp? by the equation 
V 
- 3 
V. = — (3) 
The dependence of Vz on R,,, for Vp = 6 volts and Ry = 


10° ohms, is indicated in Fig. 3b. The lowest value of Ry 
occurs when saturated liquid is flowing in the test section. 
During the course of the experiments, Ry for saturated liq- 
uid was found to vary between 10? end 2 x 109 ohms. For the 
lower 110115 Of R,, = 10° ohms, it can be seen that for Ry = 
10° ohms, the maximum sensitivity will occur in the liquid 
region. For Ry = 00 the maximum sensitivity will occur at a 
much higher Ro corresponding to a flow pattern with a nigher 


void fraction. Consequently, by adjusting R the maximum 


y? 
sensitivity can be chosen to correspond to a flow regime 
boundary. However, since these boundaries are not known in 
advance, it is desirable to maintain R. at a fixed value for 
each run. It was found most convenient to adjust Ry so that 
maximum sensitivity was obtained in the bubbly region. This 
could not always be achieved due to the restriction imposed 

by the internal resistance of the oscilloscope. This, however, 


did not affect the results to any noticeable degree since the 


variable of interest was the shape rather than the value of VT。 


2.3 Probe Interpretation 
The output of the probe for typical runs is described 


below. Throughout the interpretation phase of the experi- 
ment, Polaroid photos were simultaneously taken of the exit 


flow and of the oscilloscope trace. 
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2.3.1 Velocity Effect on Probe Signal Amplitude 

With stationary water in the test section, a large 
Signal was obtained. When the pump was started and water 
began to flow, the signal decreased by one to two volts. 
As the water velocity increased, the voltage decreased, 
although the decrease was much less than that observed at 
startup. Other experimenten >? also observed this effect 
and speculated that this is an "electrokinetic" effect, 


associated with the charge-carrying capacity of the water. 


2.3.2 Temperature Effect on Probe Signal Amplitude 

When the power was applied to the test section, the 
bulk fluid temperature increased and the voltage trace in- 
creased. This voltage increased to a maximum as the satu- 
ration temperature was reached. This was expected since 
the conductivity of the water increases with temperature, 
causing a greater voltage drop across the variable resis- 
Cor, Ry ° To insure that this effect was only due to temp- 
erature, power was shut off and the water temperature was 


increased by means of the preheaters. The same effect was 


observed. 


2.3.3 Flow Regimes Related to Probe Signal 


The number of distinct flow regimes of two phase flow 
was shown to be somewhat arbitrary. Four main flow 
patterns were recognized in the present experiments: 
bubbly, slug, annular, and wispy annular flow. The physi- 
cal or geometric description of the observed flow patterns 


is very similar to that given in Section 1.2. 
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Bubbly flow is characterized by a continuous liquid 
phase with small vapor bubbles dispersed somewhat uniformly 
throughout the tube. The characteristic bullet-shaped 
bubbles of slug flow could not be identified under the 
present high heat flux conditions. Slug flow was recog- 
nized whenever large, irregular bubbles, whose size was on 
the order of the diameter of the tube had somewhat continu- 
ous slugs of water between them. Annular flow consisted of 
a thin film of liquid, usually quite turbulent, flowing along 
the wall, while vapor and fine droplets flowed in the center 
of the tube. This has often been called spray-annular flow, 
but for the purposes of this report. annular will be used. 
The character of the core flow was inferred from the charge 
carrying spray effect which is discussed in Sec. 2.3.4. 
Wispy annular flow was characterized by a continuously in- 
creasing void fraction without the distinct characteristics 
of the above three flows. 

Since the classification of these flow regimes is 
rather qualitative, different individuals may interpret 
the flows differently. The data in Figs. 4, 5, and 6 were 
chosen so as to indicate the semi-quantitative definitions 
of the flows which were used throughout this study. 

Figure 7 presents a summary of the probe traces for the 
various flow patterns. 

Bubbly Flow In Photo I in Fig. 4, a 2,8-volt signal with 
very small fluctuations is observed, These small fluctua- 


tions are due to bubbles which can be seen in the sight 





S102 
section. Figure 7a shows the schematic of the oscillo- 
scope trace for bubbly flow. Incipient boiling could not 
be detected by the probe due to the small amount of vapor 
generation. It could not be detected with the sight sec- 
tion either, because these bubbles recondensed in the 
subcooled bulk flow before visual observation could be 
made. 
Bubbly-to-Slug Transition The transition from bubbly to 
slug (BTS) flow was first defined when the small pips 
previously seen in bubbly flow became large enough to pro- 
duce a 1.0 - 1.5-volt deflection. Photo II graphically 
shows a large bubble just leaving the heated section. 
This bubble caused a 1.3-volt deflection on the scope 
trace. This transition, schematically shown in Fig. Tb, 
was well defined and reproducible. 
Slug Flow Slug flow was considered to be present when 
the probe signal alternated between all-liquid and all- 
vapor levels. In Photo III it is noted that clean 
bubbles and slugs, of the variety found in adiabatic flow, 
do not exist. The bubble is like a vapor void without a 
definite head or tail, and the slugs have small bubbles 
dispersed throughout. At high flow rates, the slug trace 
was somewhat different, as seen in Photo IV. The frequency 
of the full deflection is increased dramatically and the 
flow appears to be very agitated. Figures 7c and 7d show 
the first full slug and slug flow, respectively. 
Slug-to-Annular Transition The transition criterion for 


the slug to annular (STA) was more subjective than the 
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other criteria. From Pile. Td It is seen that for siug flow, 


the voltage deflections average out between (V and O. 


max 
In Fig. fe it appears that the slugs are thinner and a 
larger percentage of vapor is passing the probe. The STA 
transition was defined when the voltage deflections appeared 


to rise from O volts up,rather than from (V down. An 


I^max 
example of this ilow is seen in Photo V. It is noted that 
it would be very difficult to distinguish this flow from 
annular flow, particularly at high G, using the sight sec- 
tion alone. As test-section power increases further 

(Photo VI, Fig. 7£), the voltage trace is similar to that of 
the initial description of the STA transition except that 
the peak voltages are not full deflections. This flow ís 
still the transitional flow from slug to annular. Physi- 
cally, it appears as if most of the tube consists of annular 
flow, but much of the liquid transport is still in the very 
rapidly moving slugs. 

Annular Flow With greater power, the voltage deflections 
become smaller, thus indicating that the annular flow region 
is being approached. Annular flow was defined to be that 
flow which would give an average deflection no greater than 
0.5 volts (Fig. 7g). This definition did allow larger de- 
flections to occur, but the frequency of these large deflec- 
tions was greatly reduced. At low velocity annular flow 

was well defined since the voltage was practically zero 
(Photo VII). Photo VIII is another example of annular flow. 


The increased velocity of the vapor core caused large wave- 
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other criteria. From Fig. TA it is seen that. for slug flow, 


the voltage deflections average out between (V and 0. 


mes 
In Fig. (6 it appears that the slugs are thinner and a 
larger percentage of vapor is passing the probe. The STA 
transition was defined when the voltage deflections appeared 


to rise from O volts up,rather than from (V down. An 


1/max 
example of this flow is seen in Photo V. It is noted that 
it would be very difficult to distinguish this flow from 
annular flow, particularly at high G, using the sight sec- 
tion alone. As test-section power increases further 

(Photo VI, Fig. 7f), the voltage trace is similar to that of 
the initial description of the STA transition except that 
the peak voltages are not full deflections. This flow is 
still the transitional flow from slug to annular. Physi- 
cally, 1t appears as if most of the tube consists of annular 
flow, but much of tne liquid transport is still in the very 
rapidly moving slugs. 

Annular Flow With greater power, the voltage deflections 
become smaller, thus indicating that the annular flow region 
is being approached. Annular flow was defined to be that 
flow which would give an average deflection no greater than 
0.5 volts (Fig. 78). This definition did allow larger de- 
flections to occur, but the frequency of these large deflec- 
tions was greatly reduced. At low velocity annular flow 

was well defined since the voltage was practically zero 
(Photo VII). Photo VIII is another example of annular flow. 


The increased velocity of the vapor core caused large wave- 
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like motions on the surface of the annular film. These 
waves appeared as large deflections on the oscilloscope. 
As G was increased, these waves increased in frequency and 
the film had a milky appearance, probably due to the in- 
creased entrainment and deposition of droplets. Very 
little could have been said about the interior flow, i.e. 
waves,if only photography had been used in this study. 
Wispy Annular Under certain conditions, i.e. T,, = 200°F, 

6 


lbm/nr-ft®, a varia- 


P = 20-30 psia, G =1.6 - 2.0 x 10 
tion of the above flow regimes was noted. AS power was 
increased beyond bubbly flow, a strong bubbly trace 
travelled downward from (Vi) ax to O volts. Slug flow did 
not appear and in its place was a flow which can be called 


wispy annular ie) or frotn 9), 


The benavior of tne trace 
is schematically shown in Figs. Th, Ti, and Tj and in. the 
corresponding Photos IX, X, and XI. As the void fraction 
is increased, the resistance between probe and wall 
becomes greater until annular flow is reached. It would 
be possible to describe many flow patterns in this high- 
velocity transition between bubbly and annular fiow. 
Froth flow would probably be a more appropriate descrip- 
tion of the beginning of this transition, while wispy annu- 
lar could describe the latter part of the transition. 
Since the probe indicated that the transition was quite 
homogeneous throughout, it appears that a single classi- 


fication is adequate. 
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2.3.4 Charge Carrying Spray Effect 
Photo XII again shows annular flow, but the oscillo- 
scope trace is radically different in that negative traces 
up to one volt are seen. This type of trace was not ob- 


served in Nigh pressure studies (16,8) 


and initially caused 
quite a bit of trouble in interpretation. With the probe 
circuit battery removed and the oscilloscope connected 
Girectly across the probe, a fluctuating negative voltage 
was still observed. The magnitude and frequency of this 
fluctuation increased as the flow rate was increased, even 
though no voltage source was present. A possible pickup 
effect on the probe from the voltage across the test tube 
was discounted because the generator output voltage did 
not have such large fluctuations. 

The explanation of this negative voltage effect is 
believed to be as follows. As annular flow is approached, 
small droplets of water are torn away from the liquid 
film and are transported as a spray in the core. By some 
physical effect these droplets become negatively charged 
upon detaching from the bulk flow, thus leaving the film 
positively charged. The droplets travel down the tube and 
rejoin the film and lose their charge, escape through the 
exit and lose their charge in the loop, or deposit their 
charge on the probe. This last possibility appears to be 
the source of the negative voltage observed on the scope. 
Professor Ds. P. Ke115 (18) has observed such effects witn 


aerosol-type flows and indicated that voltages up to 
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several volts could be observed if there were enough drop- 
lets. The mechanism that causes the detached particle to 
have a negative charge has not been identified. It appears 
that this effect is quite important and should receive more 
attention. Vonnegut, et a1, 09) indicate that such charged 
water droplets are produced naturally in clouds and thunder- 
storms and as spray from lakes and oceans. 

This negatively-charged droplet theory is substanti- 
ated by several observations. When the water was particularly 
pure, i.e. 3 megohm-cm and above, this effect was quite 
strong. However, after adding salt to the water (30-40 
PPM), the negative voltage was not observed. The sketches 
k and 1 “in Fig. 7 show this purity effect. 

The more conductive water allows more of the negatively 
charged particles to transfer their charge to the liquid 
film rather than depositing it on the probe. This also 
explains why the Dynatech workers, who used a similar 
probe, did not observe this effect. They were operating 
at 500 and 1000 psia and had a correspondingly higher temp- 
erature level. The conductivity of water increases greatly 
as the temperature increases; for example, the conductivity 
of pure water due to ionization is about a 1000 times 
greater at 5469F than at 779p. (20) Consequently, this con- 
ductive water allowed the negative charge to leak into the 
liquid film before hitting the probe. with this theory, it 
would be expected that tne negative voltage would increase 
with increased mass flow rate. This effect was also obser- 


ved, as stated earlier. 





ue 
Several methods were attempted to eliminate this effect 
electronically by using a higher battery voltage, and by 
using an ac high-frequency voltage source. The first 
method reduced the relative voltage of the effect but did 
not eliminate it. By means of filter circuits and high 
frequency ac, it was hoped to eliminate the spurious nega- 
tive charge and retain only the desired signals. However, 
the filters and bias circuits also removed the desired 
fluctuations. With experience it was found that proper 
interpretations of the flows could be made even though the 


negative voltage was present. 


2.3.5 Reproducibility and Reliability of Probe Output 


About a dozen runs under various conditions were 
repeated for a second time. When these runs were compared 
to the original ones, it was seen that they were almost 
identical. It was concluded that any run could be repeated 
with the expectancy of getting very similar results. 

The sight section observations were quite helpful in 
interpreting the probe traces. However, when the bulk 
temperature was moderately subcooled, the vapor condensed 
rapidly in the sight section and it was necessary to look 
into the exit of the heated section to get an accurate 
visual indication of the flow pattern. When the void frac- 
tion was large, the flow was usually so disturbed that it 
was difficult to visualize the flow pattern. In both of 
these situations the probe was able to distinguish the 


important features of the flow. 





2.4 Test Program 

After the probe calibration and interpretation runs 
were completed, the actual flow-regime investigations were 
started. The following table lists all the sets of runs 
completed in this experimental program. 

All the data are presented on flow regime maps. The 
coordinates of mass velocity (G) and mass quality (x) were 
chosen because it appeared that they were best suited in 
presenting the various effects studied. In addition, they 
are usually the most convenient for engineering applica- 
tions. Researchers in adiabatic two-phase flow frequently 
make use of volume flow rate and volume quality coordinates. 
They could not be used for this study since nonequilibrium 
effects dominated and these quantities cannot be calculated. 
Several runs are also presented in terms of heat flux and 
quality to demonstrate a particular effect. 

Due to the manner in which the loop operates, it was 
difficult to control the exit pressure very accurately. 
Consequently, pressure comparisons were made at low 
pressures where the exit pressure varied from 20-35 psia 
and at high pressures between 90-100 psia. These small 
pressure differences did not affect the results to any 
appreciable degree, in comparison to the effects due to the 


large difference between high and low pressures. 


2.5 Data Reduction 
All data were reduced by machine computation. The IBM 


7094 Computer at the MIT Computation Center was used. The 
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program is described in Appendix C. The heat flux was 


determined by the formula 
_ El 
(q/A) - TDL (4) 
The equilibrium exit quality is given by 


q/w + با‎ 5-7 


: )5( 


x 
fg 
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Chapter III 
PRESENTATION AND DISCUSSION OF RESULTS 
3.1 Effects of Variables on Flow Regimes 

Experiments were conducted to investigate the effects of 
velocity, inlet temperature, exit pressure, quality, purity, 
length, diameter and orientation on the flow regime boundar- 
ies. Figures 8 through 20 present all of the flow regime 
data in terms of mass velocity versus quality. In accordance 
with the discussion in Chapter II, the following flow patterns 
were identified: bubbly, bubbly-to-slug transition, slug, 
slug-to-annular transition, annular, and wispy annular. 

It is seen in all the maps that the bubbly and slug 
flows occurred primarily in the subcooled region, while full 
annular flow was developed at 2 to 4 percent quality. The 
two-phase flow is therefore in a highly nonequilibrium state 
and a considerable volume of vapor coexists with the sub~ 
cooled liquid. 

From all the maps, it is observed that the transitions 
occur at greater subcooling or lower quality as the velocity 
of the flow is increased. 

The influence of the remaining variables on the flow 
regime boundaries can best be illustrated by composite plots. 
At this point it is desirable to consider only the major 
flow regimes. It is noted that the transition regions 
occasionally cover a broad enough range of quality so that 
they could possibly be classified as flow regimes. However, 


practical considerations dictate that the number of regimes 
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be kept to a minimum. Boundary lines were drawn through 
the center of the transition regions and these lines were 
then transferred to the appropriate composite plots. These 
plots are presented in Figs. 21 through 27. 

Figures 21 and 22 indicate that the transitions occur 
at higher quality as the inlet temperature is increased 
from 100 to 200?F. The transition bands for T, = 200°F 


are much narrower than those for T, = 100°F (Figs. 10 and 


i 
12). In addition, it is noted that slug flow was not ob- 


6 © with the higher inlet 


served above G - 1.2 x 10" lbm/hr-ft 
temperature. These latter effects were also observed for 
the smaller diameter test section (Fig. 19), but were not 
observed at higher pressure (Fig. 1H). 


At T, = 100°F, at 30 and 100 psia, particularly at 


i 
high G, the transition from slug to annular flow was notice- 
ably more violent causing the test section, the connecting 
rubber hoses, and piping to vibrate considerably. While 
testing the L/D = 30 tube, violent vibrations were also 
observed at this transition. This tube burned out in this 
region in spite of the fact that the heat flux was not 

near the critical flux reported by Lopina . (10) 

Pressure had a considerable effect on the location of 
the flow regime boundaries. Increasing the pressure from 
30 psia to 100 psia shifted the flow regimes to substan- 
tially lower quality. Figure 23 shows this pressure effect 
at T, = 100 and 200°F. 


In order to eliminate the negative voltage problem 


previously described, salt was added to increase the conduc- 


w -- 
"E u 
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tivity of the water. The effect of this dissolved impurity 
on the flow regimes is illustrated in Fig. 24. The BTS 
transition for very pure water (2.5 megohm-cm) occurs at 
slightly higher quality than the transition for the 1 and 
9 PPM water. There was no discernible effect in the STA 
transition. Conductivity variations between 1 and 3 
megohm-cm produced no noticeable shift in the flow regimes. 

Length-to-diameter ratios of 30, 60 and 90 were 
examined. From Fig. 25, it is seen that the flow regimes 
are quite similar, although for an L/D = 30 the transitions 
are at slightly lower quality. 

For the same m there was no discernable difference in 
the flow regime maps for the two diameters (Fig. 26). 

In Fig. 27 it is seen that the transition lines have 
the same general shape for the horizontal and vertical test 
sections. There is a slight shift to lower quality for the 


vertical section. At velocities below 0.4 x 10° 2 


lbm/hr-ft 
slight stratification was visually observed for the bubbly 
and BTS transitions. Since the vertical test section 
showed only a minor flow regime variation from the horizon- 
tal section, it is evident that this stratification did not 
affect the probe output at the particular position, i.e. 
center, that the probe was in. A modification of the 
probe design, whereby the probe could be moved radially 
across the tube during operation, would allow examination 


of the flow very close to the walls and thus would detect 


any stratification. 
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Explanation of Results‏ 3.2 
The various effects noted above appear to be very‏ 
closely related to incipient boiling, heat flux, and tur-‏ 
bulence of the flow. To explain this relationship, the‏ 
requirements for the transition to slug flow are first‏ 
examined.‏ 


Radovich and Moissis ^.) 


have shown that the frequency 
of bubble collisions is very high in two-phase bubbly flow. 
A relatively small number of these collisions results in 
coalescence. As more bubbles coalesce, the bubbles 
increase in size to such an extent that slug flow is obser- 
ved. The time that it takes to cause the first identifiable 
existence of slug flow has been called the agglomeration 
time. The agglomeration time is thus determined by the 
collision frequency and the probability of coalescence per 
collision. The latter quantity was found to be strongly 
dependent on the purity of the liquid. Since the liquid 
used in this study was maintained at relatively constant 
purity for most of the experiments, the agglomeration time 
depended primarily on the collision frequency. A decrease 
in the agglomeration time would mean a shifting of tne BTS 
transition to lower quality. This could be accomplished 
by increasing the number of bubbles and by causing greater 
turbulence in the flow. 

The requirements for incipient nucleation in a forced- 
convection surface-boiling system have been considered in 


detail by Bergles and Ronsenows see? When their correlation 





= 
was applied to the present conditions (Appendix D), it was 
found that incipient boiling occurred at greater subcooling 
as G increased, T, decreased, and Pa increased. Jfurther- 
more, it can be seen from Figs. 28 to 31 that the heat flux 
level increases as G increases, m decreases,. and 1 
increases. As the heat flux increases there is an increase 
in the bubble population. Bowring (23) has concluded that 
as neat Tiuk and velocity increase, there is an increase 
in the subcooling at which detatchment of the bubbles from 
the heated surface occurs. This íis due primarily to a 
general increase in turbulent mixing under these conditions. 
This enhanced mixing should promote the agglomeration of 
bubbles. 

Thus, on the basis of incipient boiling, heat flux, 
and mixing considerations, one would expect a shift in the 
BTS transition to greater subcooling as G increases, T. 
decreases, and P. increases. These trends are observed in 
the present data, as was discussed in the preceding section. 
Since the heat flux level increases as the tube length 
decreases, the shifting of the BTS transition to greater sub- 
cooling as length is decreased is also accounted for. 

Regarding the other transitions, it is reasonable to 
expect that if the BTS transition occurs at lower quality, 
the STA transition should also occur at lower quality. In 
particular, the increased bubble agitation at the higher 


heat fluxes would be expected to break up the slugs, thereby 


promoting annular flow. 
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The limited occurrance of wispy annular flow can be 
related to the above arguments. It is recalled that this 
regime occurred only for high inlet temperature and low 
pressure. Under these conditions the heat flux level is 
low and there are fewer bubbles and also less tendency for 
the bubbles to agglomerate into slugs. The core flow 
remains bubbly or forthy in appearance until the void 
fraction is so high that the spray annular configuration 
is reached. In general the wispy annular regime is only 
observed at high velocity where the short transit times are 
less condusive to formation of slug flow. 

At this point it is worth mentioning that Radovcich 
and Moissis predicted that the BTS transition would be much 
more rapid for pure water. As discussed in the preceding 
section, the present tests indicated that this transition 
for the demineralized water occurred at about the same 
quality as that for the water with dissolved salt. 

Their observation was for adiabatic flow; thus, with heat 
addition the other mechanisms appear to overshadow the 


impurity effect. 


3.3 Interpretation of Critical Heat Flux Data by Means of 


Flow Regimes 


The present flow regime observations should be useful 
in interpreting heat transfer and pressure drop data for 
boiling water at low pressures. As an example of the 
application of the present data, one of the important 


effects in bulk boiling burnout will be considered. 
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Recent studles of bulk boiling at low pressure have 
confirmed that the critical condition coincides with the 
disappearance of the liquid fiim from the heated reee) 
The film can be depleted by evaporation and entrainment, or 
it can be disrupted by bubble nucleation. In addition to 
these local effects it appears that unstable upstream flow 
patterns can produce a disrupticn of the annular film 
further downstream. 

Lopina (20) found that the burnout heat flux increased 
as the inlet temperature increased. Typical data illustra- 
ting this effect are plotted in Fig. 33. The flow regimes 
for these conditions are superimposed on the operating lines. 

First consider operation at high inlet temperature as 
shown in Fig. 33. The majority of the tube is in annular 
flow with only moderate increases in heat flux. The slug 
flow region of the tube becomes increasingly violent, but 
also moves toward the iniet, as the heat flux is increased. 
This unstable slug flow region should have little effect on 
the critical heat flux, which for these flow conditions is 
probably caused by nucleation-induced breakup of the film. 
Now, as the inlet temperature is decreased, the siug flow 
region remains in closer proximity to the exit of the tube 
and the disturbance created within this region would be 
expected to propagate downstream and affect the film at its 
thinnest point, which is at the exit of the tube. A 


significant reduction in the critical heat flux is thus 
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observed as the inlet temperature is increased. 
with short test sections the operating region will lie 
in the vicinity of zero quality. In this case the unstable 
slug flow regime is at the exit to the test section and 
relatively low burnout fluxes, noted in the present L/D = 


30 tests, for example, would be expected. 
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Chapter IV 
CONCLUSIONS 

4.1 Conclusions 

The conclusions of this investigation can be summarized 
as follows: 

1) The electrical resistance probe is exceptionally 
well suited for flow regime studies. It is simple, inexpen- 
sive, and provides reliable and consistent results. For 
low pressure studies, where nonequilibrium effects dominate, 
the probe is more reliable than an exit sight section. 

2) The interesting phenomenon of the spray particles 
carrying a negative charge has been observed. The cause of 
this effect is not yet known, but various researchers are 
investigating the problem. Although it influenced the 
probe output under certain conditions, it did not inhibit the 
flow regime analysis. 

3) Flow regime observations were made for low pressure 
water over a wide range of flow conditions. The major flow 
regimes observed were bubbly, slug, annular, and wispy annu- 
lar flow. The flow regime maps show that bubbly, transition 
to slug, slug, and transition to annular flow occur in the 
subcooled region under nonequilibrium conditions while a 
fully developed annular flow exists by the time 2-4 percent 
quality is reached. 

4) The observed effects of the parameters on the 


flow regimes are summarized in the table below. 





I 
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Table II 
Direction of 

Change in Direction of Shift of 
Parameter Parameter Flow Regime Boundaries 
Mass Velocity Increase towards greater subcooling 
Inlet Temper- Increase towards less subcooling 

ature 
Pressure Increase towards greater subcooling 
Purity Decreased negligible change 
Length Increased slight change towards less sub- 
cooling 

Diaméter Decreased negligible change 


Orientation Horizontal to negligible change 
vertical 


5) The influence of velocity, temperature, pressure and 
length on the flow regime boundaries can be explained in terms 
of incipient boiling, heat flux, and turbulent mixing consid- 
erations. 

6) The present observations were applied to bulk boil- 
ing burnout data. It was found that the effect of inlet temp- 
erature on burnout can be related to unstable upstream flow 


patterns. 
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APPENDIX A 


Burnout Data 


19 and 21 July 1965 


= 100°F 

= 0.242 in. 

= 60 

P. G x 1076 x (a/A)op X 1076 

(psia) (1bm/hr-ft^) percent (Btu/hr-ft°) 
27.0 0,5 33:0 0.962 
28.0 0.4 37 30 0.843 
43.5 0.8 30.6 1,51 
34.0 0.8 19.0 1.12 
37,3 192 116 و1‎ 
36 .0 152 10.9 due 
51.0 1.0 10.76 1501 
45.0 1.6 8.35 1.64 
48.0 2:0 4.15 1.71 
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APPENDIX B 
Photographic Procedures 

For the photographic study, simultaneous pictures were 
taken of the oscilloscope trace and of the exit flow in the 
glass sight section. This was accomplished by shooting both 
pictures manually at the same time at the desired flow 
conditions. 

The Polaroid Oscilloscope Camera for the Hewlett-Packard 
Company Model 130B/BR Oscilloscope was used. The film was 
Polaroid ASA 3000 film and the camera settings were f 2.8 
and speed of 1/30 sec. 

The sight section was normally illuminated by a 
General Radio Co. Strobotac. The Strobotac was always off 
when actually taking photos. 

For the sight section photos, Polaroid ASA 200 Film was 
used in order to obtain better quality pictures. A General 
Radio Co. Type 1530-A Microflash was used to illuminate the 
exit test section. The 2-microsecond light pulse effectively 
stopped the motion of the flow in the sight section. Various 
backgrounds for the sight section were tried with black 
velvet giving the best results. The Polaroid camera was 
mounted on a tripod and the +2 and +4 Polaroid close-up 
lenses were used. The camera settings were f.8 and 1/30 sec. 

Below is a sketch of the lighting setup viewed from 
above. The microflash was approximately 12 inches below the 
test section and was inclined so that the light was pointed 


directly at the sight section. 
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APPENDIX C 
Computer Program for Data Reduction 

The experimental data obtained was reduced in accor- 
dance with the equations in Section 2.5 by a program written 
in FORTRAN and run on the MIT Computation Center IBM 7094 
Computer. 

The input data consisted of the run identification 
number, the exit pressure, voltage across the test section, 
shunt voltage (current through test section), enthalpy of 
the incoming fluid,and the mass flow rate. 

The program calculated the rate of heat transfer, heat 
flux, rate of heat transfer per pound mass of water, enthal- 
py of saturated liquid, latent heat, and quality. 

Below is a list of the variables as used in the program 


and their physical counterpart. 


PROGRAM VARIABLE PHYSICAL COUNTERPART 

I, J, K Numeral run identification numbers 

P EXIT Pe - the pressure at the test section 
exit (psia) 

E V - voltage across test section (volts) 

CUR I - shunt voltage (mv) 

HIN hin - enthalpy of incoming subcooled 

water 
W w - mass flow rate (lb/hr) 


A - heat transfer area (Changed this 
value for each L/D and D) (ft2 


HADDED q - rate of heat flow (Btu/hr) 


HTFLUX q/A - heat flux (Btu/hr-ft*) 





io. 


BTULBM q/w - rate of heat flow per lbm of 
water (Btu/hr-lbm) 
ENFPX hfe - enthalpy of saturated liquid at 
exit (Btu/lbm) 
ENFGPX hfg  - latent heat at exit (Btu/1bm) 
e 
X X - quality 


The enthalpy of saturated liquid (hr) as a function 
of pressure was obtained by Todreas(25) by a power series 


fit of the Keenan and Keyes' values. For P € 200 psia 


9 


1 
h.(P) = Za,P 
f is 


a. - 1.1222734(10°) 2.1239014(107!) 


© 
|| 


o 
a, = 6.3204790(10°) : - -1.0067598(107?) 
a, = -1.4742752(1077) a, = 2.9480958(1071*) 
a, = 2.5403593(107?) ag = ~4 .8430355(107*?) 
a, = -2.8788220(10°) ag = 3 .4075972(10728) 


The latent heat (heg) was obtained by fitting the 


Keenan and Keyes! values on a power series. For P <400 psia, 


2 1 
Ne, (P) = ړوغ‎ 
a, = 9.9704457(10") a, = -1.2499843(10°7) 
a, = =2.1762801(10,) ag = 1 .9322768(107*) 
a, = 1.9558791(10”°) a, = -1.8023251(1071?) 
a, = -1.2610018(10°*) ag = 9 .2680092(107*7) 
a, = 5.0308804(107! ) ag = -2 .0148643(107°*) 


The values of the heat transfer area for the differ- 


ent tubes are: 





my ee 


TUBE A (265) 
0.242 in. i.d. L/D = 30 0.038274 
0.242 60 0.076548 
0.242 90 0.114822 
0.094 60 0.0154 


The output data consisted of all the input variables 
and the computed quantities. 
The program used for the L/D = 60, D = 0.094-in. test 
tube computation follows. 
COMPUTATION PROGRAM 
1 READ 2, 1, J; K, PEXIT, E; CUR; HIN, W 


A = .0154 

HADDED = 205.36 * E * CUR 
HTFLUX = HADDED / A 
BTULBM = HADDED / W 

P = PEXIT 


5 ENFPX = 1.1222734E2 + 6.3204790*P - 1.4742752E-1*(P**2) 
1 +2.5403593E-3 *(P**3) E +2.1239014E-7*(P**5) 
2 a a +2 .9480958E-12*(P**7 
3 -4,8430355E£-15*(Px**8) Ta eae! 
6 ENFGPX=9 .9704457E2 -2.1762821*P +1.955 D E 
1 -1.2610018E-4*(P**3) +5.0308804E-7*(P**4) -1.2499843E-9*(P**5) 
2 ا‎ — a Ud 
3 49.2680092E-19*(P**8) -2.0148643E-22*( P**9 
X = (HIN + BTULBM - ENFPX)/ENFGPX 
PRINT 6, I; J KC. TEL E 
1 CUR, HIN, W, HADDED, HTFLUX, BTULBM, X, ENFPX, ENFGPX 
GO TO 1 
2 FORMAT ne 
3 FORMAT (7HORUN NOI3, 1H-, I3, 1H-,13,5X, 
OHPEXIT=E15 .6,5X, 2HE=E15..0,5X, 4HCUR=E15 .0/5X, 
UHHIN=E15 .6,5X2HW=E15.06,5X,7 
HHADDED=115 .6/5X, 7HHTFLUX=E15 .6,5X, 
THBTULBM=E15 .6,5X, 2UX=E15 .6/5X 
GHENFPX=E15 .6,5X, 7HENFGPX=E15 .6) 
END 


U FO NH 





Y 
APPENDIX D 
Incipient Boiling Analysis 

The effects of velocity, temperature, and pressure on 
the incipient boiling heat flux, (a/a), > can be examined by 
means of the correlations given in Ref. 22. This correla- 
tion for water is: 

2930 

(q/A), - 15.60 سی‎ r (D1) 

For a uniformly heated tube the heat flux can be calculated 


from: 


(q/A) 2 h(T,-7T,) 2 h (T,-T.) 4 (7,-7,) (D2) 
The saturation temperature is known from the exit pressure. 


The bulk temperature is found by the First Law: 


q = c (Tout رل‎ ل٧‎ (D3) 
which leads to: 
A 
Tout = (Age + Tin = Tb ° (D4) 


Since the probe is at the exit, only the conditions at the 
exit are of interest and therefore Ty 一 1 دا‎ The heat 
transfer coefficient, h, is found using McAdams! correla- 


(26) 


tion: 


hD, _ ap, -8,""p, 4 
Cg) - 0.023 (7) C) (D5) 
where the properties are evaluated at the bulk fluid tempera- 


ture. For the actual numerical analysis, the values of these 


properties came from Eiwakı1.(27) By assuming a value of 





- كما‎ - 
(a/a), T, and hence h can be calculated. Then using (D2), 
(T,-24,) is found. The values of (q/A) and (a Sar 


Various, GCs T ana P. can then be found and plotted. 


p 
Equation (D1) for the various pressuresis then superimposed 
on these plots. The intersections of the (q/A) vs. (T,-T.) 
curve and the (a/a), curve is the point at which incipient 
boiling is predicted. 


This analysis was carried out for the conditions that 


were presented in Figs. 28 through 31, i.e. 


6 


Gx 10 Ti E 
2 0 e 
lbm/hr-ft F (psia) 
a) 0.4 100, 150, 200 30 
b) 0.4 100, 200 100 
c) eu LOO; 150, 200 30 
d) 2.0 100, 200 100 


In Fig. 32a is a representative plot for the conditions in 
(c). To visually relate the predicted (q/A), to (x), opera- 
ting lines obtained from the formula: 


2 q/w + hin - he 


= (D6) 


£ 
fg 


are plotted on coordinates of (q/A) vs x. The (q/A), found 
from the analysis are plotted on these operating lines and are 
connected by a smooth line. This line then is the locus of 
predicted (a/a), for various inlet temperatures. This plot 
for conditions (c) above is presented in Fig. 32b. Examin- 
ing Fig. 32b it is seen that incipient boiling occurs at 


greater subcoolings for lower Ti: From the analysis of 





se 


los 
conditions a, D, d,it was also found that incipient boiling 
occurs at greater subcooling for larger pressures and, keeping 


P and T. Constant,-at higher G. 


d 
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FIG. 3a PROBE CIRCUIT DIAGRAM 
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FIG. 3b CHARACTERISTICS OF PROBE CIRCUIT 
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